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ABSTRACT 

We derive the CMB radiative transfer equation in the form of a muhipole hierarchy in 
the nearly-Friedmann-Robertson- Walker hmit of homogeneous, but anisotropic, uni- 
verses classified via their Bianchi type. Compared with previous calculations, this 
allows a more sophisticated treatment of recombination, produces predictions for the 
polarization of the radiation, and allows for reionization. Our derivation is independent 
of any assumptions about the dynamical behaviour of the field equations, except that 
it requires anisotropics to be small back to recombination; this is already demanded 
by observations. 

We calculate the polarization signal in the Bianchi VII;i case, with the parameters 
recently advocated to mimic the several large-angle anomalous features observed in the 
CMB. We find that the peak polarization signal is 1.2 fiK for the best-fit model to the 
temperature anisotropics, and is mostly confined to multipolcs I < 10. Remarkably, the 
predicted large-angle EE and TE power spectra in the Bianchi model are consistent 
with WMAP observations that are usually interpreted as evidence of early reionization. 
However, the power in i?-mode polarization is predicted to be similar to the E-mode 
power and parity-violating correlations arc also predicted by the model; the WMAP 
non-detection of either of these signals casts further strong doubts on the veracity of 
attempts to explain the large-angle anomalies with global anisotropy. On the other 
hand, given that there exist further dynamical degrees of freedom in the VII;i universes 
that are yet to be compared with CMB observations, we cannot at this time definitively 
reject the anisotropy explanation. 
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1 INTRODUCTION ing anisotropies llWaldlll983l ). this analysis is not complete 

, r , , r . , n JColiath fc Elllsl [l999l ) ■ It is therefore worth investigating 

there are a number oi observed features m the large- i i- i j- ^- i ^ i i i -i j_ 

, . . . , . . further the observationaf predictions of such modefs; whilst 

angle temperature anisotropics of the cosmic microwave , ,i i .. i., i . , r 

V J 1, 1 ^"^^ skeptic they seem unfmely, the importance of a posi- 



tive detection would be great. 



background (CMB) that are anomalous under the usual 
assu mption of statistically-isotropic, Gaussian fluctuations 

( see ICopi et al. 2007 for a recent summary). Recent analy- Us ing th e CMB models developed inlCollins fc Hawkind 

ses have suggested that a small global anisotropy in the form l|l973al ) and iBarrow. Juszkiewicz fc Sonodal l|l985h . based 

descr ibed by Bianchi models could be to blame (|jaffe et al.1 on a specific soluti on for a type-VII /i universe with only 

I2OO5I ). In classical general relativity, the dimension of ini- pressureless matter. Uaffe et all l|2005h showed that known 

tial state space leading to models with such homogeneous anomalies such as the low quadrupole amplitude, alignment 

anisotropies is always greater than the limitation of these of low-/ modes, large-scale power asymmetry and the 'cold 

models to the exactly isotropic case; furthermore, nearly spot' in t he CMB could b e mimicked. The work was ex- 

isotropic models do not nece ssarily tend to late-time isotropy tended in ljaffe et ahl (|2006l ) to include the dynamical effect 

ICollins fc Hawkind ri973bl ). Although later work has sug- of dark energy, yielding a degeneracy in the - plane 

gested an inflationary epoch can be responsible for remov- (which effectively arises through the angular diameter dis- 
tance relation, since much of the temperature anisotropy is 
generated at high redshift). The degenerate range of param- 

* Email: apontzen@ast.cam.ac.uk eters able to explain the large-angle anomalies was shown to 

t Email: a.d.challinor@ast.cam.ac.uk be inconsistent with the cosmological parameters required to 
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explain constraints such as supernovae observations and the 
CMB spectrum on smaller scales (i.e. the structure of the 
acoustic peak s). A more complete statistical analysis was 
performed by [Bridges et all (|2007l ). in which the inconsis- 
tency of inferred parameters was confirmed. 

There are remaining dynamical freedoms in the Vllh 
model w hich were not exp lored in the above papers. As 
noted in Ijaffe et"aD l)2006l ). these will need investigation 
before any firm conclusions about the compatibility of the 
model with observations can be made. 

A second important issue, hitherto not considered, is 
whether the predicted CMB polarization pattern in the 
Bianchi model that best fits the large-angle temperature 
anisotropics is consistent with current observations. The 
reason for its neglect seems to be that there are cur- 
rent ly no full predict ions for polarization in Bianchi mod- 
els (jjaffe et al.l |2006| ). In this paper, we address this is- 
sue by providing a complete and computationally conve- 
nient framework for calculating polarization in these mod- 
els, and make a first attempt at confronting the predictions 
with the three-year data f rom the Wilkins on Microwave 
Anisotropy Probe { WMAP: IPage et al^l2007^ . Polarization 
of the CMB in anisotropic cosmologies was first discussed 
by iReed l| 19681 ). where a calculation for Bianchi I mod- 
els was outline d. Further considera t ion h as been given to 
the problem in iMatzner fc TolmanI (Il982l ). where Bianchi 
types V and IX were considered. There are also a num- 
ber of papers considering the e ffect of homogeneous mag- 
netic fields on polarization (e.g. iMilaneschi fc Fabbrilll985l : 
iFabbri fc Tamburrand Il987l ). although we emphasize that 
the resulting Faraday rotation of the polarization is quite 
distinct from the purely gravitational effect under consider- 
ation here. 

This paper is organized as follows. We briefly review 
Bianchi models in Section [S] including a more intuitive 
derivation of their FRW limit than has previously been pub- 
lished. In Section [3] we derive the zero order evolution of 
the photon polarization direction along geodesies, and in- 
corporate these results into a multipole treatment of the 
Boltzmann equation that describes the radiative transfer in 
Section [4] Using this, we specialize our model to calculate 
the temperature and polarization anisotr opies expected in 
a VIIj, model with favoured parameters (|jaffe et al.ll2006l ) 
in Section [5] and discuss the difhculties of reconciling the 
results with existing observations. We defer a full statistical 
reanalysis of Bianchi signatures, including reionization, to 
future work, where we will also explore more fully the extra 
dynamical degrees of freedom in the models. 



2 BIANCHI MODELS 

In this section we give a brief review of the framework 
that we use for our calculation. Anisotropic, homogeneous 
cosmological models can be classified according to the com- 
mutation relations of their spatial symmetry groups. The 
most popular classification method is based on Bianchi's 
(1897) classification of three-parameter Lie groups. For the 
development of t h ese c la ssifications in the cosm o logica l 
context see iTaubl lll95ll): iHeckmann fc Schuckind lll962l): 
lEstabrook. Wahlguist fc Behij l|l968h : lEUis fc MacCallumI 



J 19691) ■ and for 



i^iews see IWainwrightl (|l997l ) 

and lEUis fc van Elsd (|l998l ). 

We adopt a — I — h- 1- metric signature, and will use Greek 
spacetime indices (0 — > 3) for tensor components in a general 
basis, early Latin indices (a, b etc. running from ^ 3) to 
label the vectors of any group-invariant tetrad, and middle 
Latin indices {i, j etc. running from 1 — > 3) to label spa- 
tial tetrad vectors in expressions only involving the spatial 
vectors. For the latter, we also use upper-case early Latin 
{A, B etc.) when the tetrad is time-invariant (see below). As 
usual, round brackets denote symmetrisation on the enclosed 
indices and square brackets denote anti-symmetrisation. 

A spacetime is said to be homogeneous if it can be 
foliated into space sections each admitting at least three 
linearly- independent Killing vector fields (KVFs) {^} (so 
that C(^% — 0, where £, denotes the Lie derivative and g 
is the metric tensor) with at least one subgroup that acts 
simply transitively in the space sections. We denote the ele- 
ments of the three-dimensional subgroup by where i — 1, 
2, 3. The commutator of two KVFs is also Killing and, since 
the form a subgroup, we must have 

[i^,i,]-cU,- (1) 

The C*j- = Cfjj.] are the structure constants of the Lie alge- 
bra of the group and are constant in space. (Through the 
foliation construction below they may also be shown to be 
constant throughout time.) The Jacobi identity, 

+ + li„ = 0, (2) 

restricts the structure constants to satisfy 

C^jC'ifc] ~ ^- (^) 

Since a constant linear combination of KVFs will also be a 
KVF, one is permitted to perform global linear transforma- 
tions: 



Si Si -^iSj 



(4) 



under which the structure constants transform as a (mixed) 
3-tensor. Classifying all homogeneous spacetimes amounts 
to finding solutions to equation ^ that are inequivalent 
under the linear transformations (|4]). 

The fiducial classification is achieved by decomposing 
the structure constants into irreducible vector and pseudo- 
tensor parts: 



^77 



km I rT7 

eijkn + Oidj 



n Oi — 



(Jacobi identities). 



(5) 



where n'-' is symmetric and eijk is the alternating tensor. 
Linear transformations can be used to diagonalise n*-* and, 
since Oi is an eigenvector of n*-* one may take, without loss 
of generality: 



= diag(ni,n2,n3) 



(a, 0,0). 



(6) 



By rescaling (possibly reversing their directions) and rela- 
belling the KVFs, we may reduce the structure constants 
to one of 10 distinct canonical forms t hat describe the 10 
possi ble group types (see, for example, lEllis fc MacCallumI 
Il969l ). In these forms, all non-zero m and a are either ±1 
except for the two types with an2ns 7^ in which case there 
is an additional parameter h = a"^ / {n2n-i) . 
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In order to express homogeneous tensor fields in an in- 
variant way within this space, one must construct a suitable 
basis tetrad {ea} that is invariant under the group action, 
i.e. satisfying 



(7) 



Any homogeneous tensor (T where C{T — 0) has compo- 
nents relative to this tetrad that are constant throughout 
space. The unit normal n to the hypersurfaces of homo- 
geneity (satisfying n • = 0) is necessarily group-invariant 
([^i^J = 0) viS' l-lie Leibniz property of jC, so we can al- 
ways take eo = n. Because of the foliated construction of 
the spacetime, one may label the hypersurfaces with a co- 
ordinate time t such that tIq = —t^a and hence eo can be 
shown to be geodesic. In any homogeneous hypersurface, we 
can always construct the spatial part of the group-invariant 
tetrad by making an arbitrary choice for Si at a point, and 
then dragging this fram e out across the hypersur face with 
the (see, for example, lEllis fc MacCallumll 19691 ). 

In general the elements of the tetrad will not commute 
so one has 



[ea, et,] = 7ai,ec 



(8) 



where the y^t = ifab] constant in the hypersurfaces of ho- 
mogeneity by the Jacobi identities. The rotation coefficients 
encode the covariant derivatives of the tetrad: 



(9) 



where e" are the components of the tetrad vectors in an 
arbitrary basis. The rotation coefficients are related to the 
7a6 by 



Tabc = - { — dagtc + ObQca + dcQab 
+ 7tibc + Icab — fbca) , 



(10) 



where gab = g(ea, Sb) are the tetrad components of the met- 
ric, jabc = gad'Jbc ^ud da IS the Ordinary derivative along 
the direction of e^. By taking eo = n, we have goo ~ —1, 
goi = and, since n is normalized, geodesic and irrota- 
tional, ■jai, ~ 0. The group-invariance of the tetrad restricts 
the form of the spatial components ■yfj: they are related to 
the group structure constants by a linear transformation and 
can be classified in the sa me way (and, necessarily, fall into 
the same Bianchi type: see iMacCalluml 1973). 

It is always possibl e to construct the tetrad so that 
it is orthonormal; see lEllis fc MacCaUumI l|l969l ) for de- 
tails. However, a convenient alternat i ve is th e time-invariant 
frame used by ICoUins fc Hawkind (|l973ah . In this case, 
the tetrad is constructed so that [n, e^] = 0. The Jacobi 
identities applied to n, Sa and show that this is still 
consistent with the tetrad being group invariant. A spe- 
cific construction is to start with a group-invariant tetrad 
over a hypersurface and to drag the tetrad along the nor- 
mal n to cover spacetime. An important consequence is 
that the are constant throughout spacetime for the 
group-invariant tetrad. Moreover, 7^^ = and the only 
non- vanishing components are spatial, i.e. 7,* . Because the 
Bianchi classification of the Killing group and tetrad com- 
mutators is the same, it is always possible to perform a 
global (three-dimensional) linear transformation of the time- 
invariant tetrad to bring the equal to the canonical struc- 
ture constants, i.e. 7* = Cij. We shall use upper-case early 



Latin indices {A, B etc.) to denote the spatial elements of 
this specific time-invariant tetrad. 

The spacetime metric can then be written in terms of 
the time-invariant tetrad as 



(11) 



36 n satisfy e^eg 



g = — n ® n + gAB (i)e' ® e^ 
where the basis one-forms 

5g. Here, as usual, g^^ denotes the inverse of gAB- The 
e'* are also group- and time-invar iant. We deco mpose the 
spatial part of the metric following iMisnej l|l968l) : 



2a(t) I 2i3(t)\ 

gAB = e ^ ' [e^^ ') 



(12) 



where /3 is a matrix with zero trace and the matrix expo- 
nential is defined as usual via Taylor expansion, yielding 
dete'' = 1. Thus a represents shape-preserving expansion 
whilst /3 represents volume-preserving shape deformation. 
The expansion rate of the n congruence is V^n" = 3(i and 
its shear - i.e. the trace-free, symmetric spatial projection 
of V^n^ - has components 



CAS 



(13) 



Here, and throughout, overdots denote derivatives with re- 
spect to t. We shall make use of conformal time, rj, defined 
by dt — e"dr], and we denote derivatives with respect to 
with primes. 

The field equations are generally more naturally ex- 
pressed in a (group-invariant) orthonormal frame, since the 
metric derivatives then vanish and one obtains first order 
equations for their commutation functions. The spatial vec- 
tors of such a frame will be represented with lower-case mid- 
dle Latin indices, i, j etc. One can al ways define a specific 
orthonormal frame by (|Hawkinell 19691 ) 



e-A- 



(14) 



The orthonormal tetrad is completed with n. The compo- 
nents of the shear in the orthonormal frame are 



(e'') fe(i(e '^)i)k 



(15) 



Note that this corresponds to the matter shear tensor only 
in the case that the fluid flow is not tilted relative to the 
hypersurfaces; i.e. the case n = v where v is the fluid 4- 
velocity. In this paper, we do not specialize in this way. 



2.1 FRW limit 

It is instructive to sketch the derivation of the FRW limit 
of spacetimes in the time-invariant fram e (the calculation in 
the orthonormal frame is described in lEUis fc MacCaUumI 
1969 and yields identical conditions). 

A homogeneous spacetime can only be FRW if the ex- 
pansion and 3-curvature are isotropic. Vanishing shear re- 
quires /3as = const, and then we can always set /3as = 
(i.e. diagonalise the spatial metric) with a constant linear 
transformation of the &a- After the transformation, the 7*^ 
will no longer be equal to the canonical group structure con- 
stants (so we denote them with lower-case indices) but we 
can perform a further (constant) orthogonal transformation 
of the ei (hence preserving /3 = 0) to bring the to the 
canonical form up to positive scalings of a, ni, 712 and 713. 
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Type 


a 


ni 


r»2 






I 

















V 


-1 











-6 


VIIo 








1 


1 





VII,, 







1 


1 


~eh 


IX 





1 


1 


1 


3/2 



K = E{n+p'e^) 



(20) 



Table 1. Bianchi groups with FRW limit and their structure con- 
stants in canonical form. In the case Vllfc, we take a = — \/fen2n3 
for compatibility with th e notation of lCoUins &c HawkineJ lll973al ) 
andU arrow et al. I lll985l) . Note that a can also be positive, but 
this case will be related by a rotation of the sky (constructed 
by ei — ► — ei followed by 63 — > —63) and need not be consid- 
ered separately. However, the overall parity must be considered; 
one way is to consider the transformation ei — > — ei which re- 
verses the sign of a and all rii . The final column is the comoving 
3-curvature scalar in the /3 = FRW limit. 



With P = 0, the 3-curvature is 



(3) 



-■ylnli + '2'ykiiYji + la)] 



(16) 



Zero- and first-order expressions (such as the one above) will 
often contain surprising index placement; note that metric 
factors are always included where necessary to account for 
this. 

In terms of the decomposition of 7^* , one obtains the off- 
diagonal element '■^'i?23 = o.{n2 — W3), yielding the condition 



a = or 712 = ns 



(17) 



if the curvature is to be isotropic. The vanishing of the re- 
mainder of the trace-free part further requires 



nf + 712723 = 7I2 -I- ni7l3 = + 711712, 



(18) 



which implies that at least one of the Tii are zero and the 
remaining two are equal. By examining a complete list of 
distinct group types, these conditions yield the spaces with 
an FRW limit (Table [T}. In each case we can use a volume 
rescaling (change of q) to set the non-zero rii to unity or, 
if the rii vanish, a to —1. This is enough freedom to set 
the 7,* of the time-invariant frame, with 13 — 0, equal to 
the corresponding canonical group structure constants. For 
a universe that is close to FRW, we can therefore treat (3 
as a small perturbation while working in the time-invariant 
frame with 7^^ still in canonical form. 



3 PHOTON PROPAGATION 

We parameterize the photon pro pagation dire ct ion v ector 
according to the convention in iBarrow et all l| 19851 ); al- 
though this unusually places the azimuthal direction along 
ei, it more conveniently reflects the symmetries of the fldu- 
cial Bianchi classification, which ultimately yields much sim- 
pler equations. Thus in the orthonormal frame, the photon 
direction has components 



p — (cos 0, sm a cos cp, sm sm cp) 
and the photon 4-momentum is 



(19) 



where E is the photon energy measured by an observer on a 
hypersurface-orthogonal path. It is convenient to introduce 
the comoving energy, e = Be" . In terms of this, the spatial 
components of K on the time-invariant tetrad (in which the 
geodesic equations take their simplest form) are 

KA-e{e^)^^p. (21) 

The energy change along a geodesic follows from differ- 
entiating E = K ■ n. The geodesic equation then gives the 
exact result 



-ee p p'Gij, 



(22) 



where e = de/drj. This is e quivalent to the integral result 
obtained by iHawkind l|l969l ). For the evolution of 6 and 0, 
we make use of the (exact) geodesi c equation in the time- 
invariant frame (|Barrow et al.|[l985l ): 



KcKd „c 



(23) 



We only require the evolution of 6 and (j) to zero order in /3 
since the radiation is necessarily isotropic in the FRW limit. 
Setting /? = in equation (|23|l . we find 



[a + (713 — 712) cos 4> sin (/>] sin 6 

[711 — 713 + (W3 — 7I2) COS^ (j)] COS 6. 



(24) 
(25) 



We denote the (complex) polarization 4-vector by P. 
It is normalized so that P ■ P* — 1 and, in the Lorentz 
gauge, is orthogonal to K. The polarization P is parallel 
transported along the photon geodesic. We are more inter- 
ested in the observed polarization P (i.e. the electric field 
direction for radiation in a pure state) relative to n. This 
is given by screen-projecting P perpendicular to n and the 
photon direction p: 

P"=H"pP^, (26) 
where the screen-projection tensor TLap is defined by 

HaH = gal3 + naUp — PaPl3- (27) 

The e volution of P fo llows from the parallel transport of P 
itself (|Challinoj|2000h : 



0. 



(28) 



We find for the tetrad components of P in the time-invariant 
basis the exact equation 

Pa - P^'aAB -PaP^'p'^obc - dpA + p" P'^Tabc = 0, (29) 

where Fabc are the rotation coefficients given by the spatial 
components of equation (|10[) FI There is only one degree of 
freedom in the projected polarization vector, described by 
the angle between it and the ee direction; 



P = eg COS tp + e^ sin t/j 



(30) 



where eg and are constructed from the orthonormal- 
frame vectors. The components on the orthonormal frame 
are therefore 



^ The metric derivatives vanish in this case since gjj is constant 
in each hypersurface. 
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C— sin 9 cos 
cos 9 cos (f> cos ip — sin sin tp 
cos 9 sin (ji cos ip + cos sin tp 



(31) 



The polarization is first-order in /3 so we only require the 
evolution of ip at zero order. In this FRW limit, equation (|29[) 
reduces to 



2tp' — —Til — (ns — 712) cos 20. 



(32) 



Since we are considering models close to FRW, 713 = 712 
(see Table [ij and our evolution equations simplify to 

9' = -'Vhsm9 
4>' = (ni — 713) cos6' 
2ip' = -711, (33) 

where, for type V universes, one takes ^/h — 1. We see that 
the polarization only rotates relative to eg and in type-IX 
universes. These general equations for 9 and (j) agree with th e 
specific cases given in Appendix B of iBarrow et al. 
the general tp equation ag rees with the previous ana lysis of 
type-IX universes given in iMatzner fc TolmanI (|l982l '). 



4 BOLTZMANN EQUATION 

In this section, we derive the first-order Boltzmann equation 
describing polarized radiative transfer in any Bianchi model 
with an FRW limit. 

CMB polarization is generated by Thomson scattering 
and is therefore expected to be only linearly polarized. We 
can describe the radiation distribution function in terms of 
Stokes parameters /, q and u where / gives the expected 
number density of photons per proper phase-space volume 
irrespective of their polarization state. The parameters q and 
u describe linear polarization relative to a basis that we take 
to be ee and for propagation along p\ Then (/ -f- g)/2 
is the expected phase-space number density of photons that 
would be found in the eg linear-polarization state, and (/ -f 
u)/2 is the same for the state rotated by 45deg (towards 
e,^). In the absence of scattering, / is conserved along the 
photon path in phase space and q and u are conserved if 
referred to bases that rotate like P in equation (|28|l . We 
parameterize /, q and u by the photon-direction angles 9 
and (f> and by comoving energy e. Since these parameters are 
defined relative to a group- invariant tetrad, / is independent 
of position in the hypersurfaces of homogeneity. The same 
is true for the Stokes parameters since the basis on which 
they are defined is constructed in an invariant manner. 

The Lagrangian derivative of / in phase space is 



D/ 9f df ,df , df , 



(34) 



and is only non-zero because of scattering. The quanti- 
ties e', 9' and (j)' ''■re given by equations (|22|l and (I33|) . 
Anisotropies are formed through the last term in equa- 
tion (|34|l and by Thomson scattering off electrons with a 
non-zero peculiar velocity; both effects are first order. The 
remaining terms 'advect' the resulting pattern on the sphere 
using the zero-order transport equations. The energy de- 
pendence of the anisotropies is therefore proportional to 
edf/de where / is the Planck distribution function in the 
FRW background. We define the dimensionless temperature 
anisotropies Q{9, 4>; rf) by 



(35) 



Polarization is generated by scattering the anisotropies and, 
since Thomson scattering is achromatic, the polarization has 
the same e spectrum as the anisotropies. We can therefore 
introduce dimensionless 'thermodynamic-equivalent' Stokes 
parameters Q{9,(f);rf) and U{9,(f);rj) as 



{q±iu){K-n) = 



d/ 
dine 



{Q±iU). 



(36) 



The quantities Q±iU have spin-weight 2, i.e. under a change 
of basis 



Q±iU 



(0±i[/).(37) 



Expressed in this way, the time dependence of the tem- 
perature anisotropies obeys: 



(38) 



where the first term on the right describes the Thomson 
scattering kernel (Section I4.3|l . the next two describe ad- 
vection of the patterns on the sphere, and the final term is 
gravitational redshifting in the anisotropic expansion due to 
the shear. The advection terms are more transparent when 
interpreted as being due to the spatial dependence (relative 
to a parallel-propagated rather than time-invariant basis) 
transforming to angular dependence through free streaming. 
For polarization, we have 

d{Q±iU) ^ D(Q±it/) _ d{Q±iU) , _ d{Q±iU) , 
dr) ^ D7? 89 d(t> 

±2i{Q±iU)tp', (39) 

where the last term arises from polarization rotation. Again, 
the first term on the right describes Thomson scattering (see 
Section 143)) . 

We expand the temperature anisotropies in terms of 
spherical harmonics about the propagation direction: 



e(p) = ^ervr(p) 

Im 



(40) 



For polarization, we expand in spin-weighted spherical har- 
monicfl as 



{Q±iU){p) 



Im 



{Er±iBr)±2YT{p). 



(41) 



Note that the Bianchi models are not parity invariant; the 
mirror universe can be obtained using the standard trans- 
formations 



er 



(-i)'er 

(-i)'sr 



(42) 



or, equivalently, fiipping the signs of a and all rn (which 
inverts the direction of the ei axis). 



^ For a brief review, see Appendices A-C 
of lLewis. Challinor fc Turokl ll2002h . 
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4.1 Gravitational redshifting 

The effect of the shear on the observed temperature pattern 
reads 



96 



-e pp'aij, 



(43) 



which is a pure quadrupole. In terms of spherical harmonics, 



del 

drj 

del 

drj 

del 

drj 
drj 



(0-33 — (722 + 2io-23) e° 
((T12 — icris) e" 



= (-ir 



-CTiie 



905 



drj 



(44) 



where aij are the components of the shear in the orthonor- 
mal frame. We see that shear injects power at Z = 2, but 
subsequent advection generally transports this to higher I. 

4.2 Advection equations 

The advection part of the Boltzmann equation for the tem- 
perature is 



I dv 



de{e,<t>)^, de{e, 



de 



(45) 



with 6' and 0' given by equations (|33p . Making use of the 
results 



Z + 1 



..a.yr I / [(/ + l)^-m2][(/ + l)^-.^] 

Smt? — — -\l . ^ s-Ii + l 



(2/ + l)(2; + 3) 



i{i + iy 



ms l + l HP - m?){P - s'^) 



COS 6*3 V;" 



I V (2/-l)(2/ + l) 
1 /[(/ + l)2-m2][(Z + l)2-s2] 



(2Z + l)(2; + 3) 



^I(7TT)^ ' +TV (2/-i)(2; + i) ^^'-1 ^^^^ 

which follow from standard recursion relations 
for the related Wigner functions - DL^s( (/), 9, 0) 



(e.g. IVarshalovich. Moskalev fc Khersonskiil Il998r ). we 
find for the advective contribution to the time derivative of 
the multipoles 

i+i 



V drj 
Here, 

/m 
1,1-1 

/m 
1,1 

/m 
1,1 + 1 



advcc i' = i_l 



P — m2 



(2Z-1)(2Z + 1) 



(; + l)2-m2 
(2; + l)(2Z + 3) 



[imAn + (/ - l)Vh\ 



imAn — {l + 2)Vh\ (48) 



with An = 713 — JT-i- Note that ±fiii^ = ^fui* , as required 
by the reality of e{d,(p) in equation (|47|l . 



For polarization, we have 
fd{Q±iU)\ _ d{Q±iU)^, d{Q±iU),, 

[ — d^ — } = de — ^ d^'^ 

\ ' / advcc 

±2i'iP'{Q±iU). (49) 

We write 
fdiEr±^Bn 



V dr/ 
which may be expressed 

dET 

drj 



dBT 

drj 



i+i 

- E 

advcc = 



±gir'(i5r ±iSr), (50) 



2 E [(+<?«' + 


m 

-Qw 


)£,? 


+n+5'ii' - - 


m - 

-9w. 




1 

^ E [(+^?"' + 


m 


)s" 


■ 1 771 

-^{+9ii' - - 


-9ii' . 





(51) 

We see that there are two modes of propagation; one trans- 
fers power amongst I and the other mixes E and _B-modes. 
The mixing terms are 



+91,1 



= — 2mi 



4 {^/h — irnAn^ m 



l{l + l) 

mm m m r\ /r(~»\ 

+9i,i+i- -9i,i+i = +9i,i-i - -9i,i-i = 0. (52) 

Polarization is generated from Thomson scattering as a pure 
electric quadrupole (see Section|33J but B-modes can subse- 
quently be produced through advection. This happens in all 
nearly-FRW Bianchi models except type I. The ni term in 
equations (|52p arises from the polarization rotation tp'; the 
remaining terms are from the evolution of the photon direc- 
tion relative to the invariant frame. For the power transfer, 
we find 

+9i,i + -9i,i = 0, 

m I m 
+Ql,l-1 ~r -9l,l-l 



" P(2l-l){2l + l) 

[{I - l)Vh + imAn] , 

™ , ™ _ „ / [(/ + l)2-m2][(; + l)2-4] 
+91,1+1 + -9 1,1+1 - (/ + l)2(a + l)(2/ + 3) 

[-{l + 2)Vh + imAn\ , (53) 

which are not affected by polarization rotation. Note that 

±9w = ^9iv , (54) 

as required by the reality of E{9,(f>) and B{6,(f>) in equa- 
tion (|50)) . 



4.3 Scattering equations 

We use a standard Tho mson scattering kernel in the form 
deriv ed bv lHu fc White! (il997i) (see also iDautcourt fc Rose! 

[lizi): 

D(Er ±iBr)\ ,f , m 

J = r'{-{Er±iBr) 
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= r' -er(i-5« 



(55) 



where r' = 
time, and 

u-i = 
Mo = 
ui = 



n^ate" gives the scattering rate in conformal 




(56) 



are the dipole moments of the electron peculiar velocity in 
the orthonormal frame. It follows from equation (I55|l that 
Thomson scattering of the temperature quadrupole gener- 
ates polarization that is an i5-mode quadrupole. 



5 IMPLEMENTATION FOR VIIh UNIVERSES 
5.1 Field equations 

The complete set of field equations are available from, for ex- 
ample, [Wainwright (1997). Naturally, these reduce at zeroth- 
order to the standard Friedmann and acceleration equations, 
so that 
a drj -1 
dz 



H{z) 



-1/2 



(57) 



where ^m.o, ^^a,o and n^.o have their usual meanings, 
(1 -|- z)~^ = e°'~°'° and H = a. The Friedmann constraint 
equation relates the group parameter h to the curvature (see 
Table [Ij: h = Hie'^°">nK,o- 

The evolution of the shear is required at first order and 
is provided by the trace-free part of the spatial evolution 
equations. 



7ij — 3H(7ij 



(3) 



in the orthonormal frame. Here, 



(3) 



Si' 



(3) 



(58) 



(59) 



is the trace-free part of the intrinsic 3-curvature of the ho- 
mogeneous hypersurfaces. Also, we have assumed a perfect 
fiuid (i.e. zero anisotropic stress) and the shear equation (|58p 
holds only at first order. It is convenient that, in the VIIh 
case, and Svi are zero to this accuracy, and, further- 

more, no coefficient of /3i2 or /3i3 enters into the expression 
for '^'iSij, so that one may study a simple model in which 
aij = except for (Ti2,(Ji3 oc e~"^". The linear constraint 
equations show that the matter in such a model contains 
vorticity, i.e. the separation of neighbouring particles rotates 
relative to inertial gyroscopes. 

Whilst it is not prohibitively difficult to implement a 
numerical solution for the most general case, we defer such 
a treatment to a later paper. Instead, we take advantage 



of the simplified so lutions to derive the polarization in the 
favoured models of Ijaffe et all l|2006l ). 

We shall adopt the standard assumption that the CMB 
signal from global anisotropy adds linearly to that from inho- 
mogeneities. This is clearly correct insofar as the linear-order 
perturbations are concerned; however, given that generic 
anisotropy modes grow towards the initial singularity, there 
is no guarantee that standard inflationary mechanisms for 
generating inhomogeneities can be invoked. Ignoring this po- 
tential inconsistency is pragmatic, but investigation would 
certainly be necessary if the resulting models gain any sig- 
nificant observational support. 



5.2 Tilt decay 

To be consistent with our assumption of small departures 
from FRW symmetry, we assume that all peculiar velocities 
are small. If we write the total momentum density of all 
matter and radiation as p('°'' = ^ (P(n) + P{n))u(n), the 
linear constraint equation relates the spatial components to 
the shear: 



8nP, 



(tot) 



(60) 



to first order in the orthonormal frame. Here, the Cj^, are 
the structure constants in canonical form. In the restricted 
solution (712, (Ti3 oc e~^°' , so, assuming the total momentum 
density is dominated by a barotropic fluid with equation of 
state p — wp for constant w, we have p oc e"'^'^"'"™^" and 

\u\ — [uiU ) ' oc ' 

constant lu = 1/3 
e"" w = G. 



(61) 



This behaviour of |u| is consistent with momentum con- 
servation. To see this, consider the Euler equation for a non- 
interacting ideal fluid in the time-invariant frame. The fiuid 
pressure is constant on surfaces of homogeneity but gra- 
dients proportional to puA appear in the fluid rest frame. 
These accelerate the fiuid so that 



puA + (p + p)uA 



0. 



(62) 



Solving gives ua oc e and, recalling the zero-order metric 
g^^ = e~^°, we recover equation (|61|l . 

The above introduces a complication in multi-fiuid 
models, which appears to have been overlooked in recent 
work. For two components, say, the tilt velocities W(i),U(2) 
need not be the same (except in the case of strong cou- 
pling). Only the total momentum density, (pi -|-pi)w(i) -f 
{p2 +P2)it(2) is constrained by the shear so there is addi- 
tional freedom in the solution|j Although the dark matter 
density will be dominant around the time of recombination, 
z ~ «LSS, the tilt velocity of relevance for the CMB is mani- 
festly that of the baryons. Given that the baryons are tightly 
coupled to the photons until the last scattering surface, they 
experience a significant pressure and their tilt decay will be 
halted; this will not be the case for the dark matter. Thus, 



We note that, with this effect in mind, the application of the 
term 'universal vorticity' to describe VII^ cosmologies is an over- 
simplification, since the vorticity of the dark matter need bear no 
resemblance to that of the baryons. 
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one needs to consider with care how to estimate the electron 
velocity in the Thomson scattering terms ((56}. 

For most of cosmic history before recombination, the 
baryon-photon plasma has an equation of state parameter 
close to w = 1/3. A simple approximation is thus obtained 
by assuming the baryon tilt velocity remains constant be- 
fore recombination, after which it decays as the inverse scale 
factor. If dark matter decouples at 2dm, the ratio of baryon 
to dark-matter peculiar velocities at last scattering will be 



t(b) 



|W{c)| 



1 + ZPM 
1 



2LSS 



(63) 



We take dark-matter dec oupling to be at redshift (see e.g. 
iLoeb fc Zaldarriaga|[2005l ) 



2dm 



lOMeV 



VlOOGeVy VlOOGeVy 



(64) 



where ks is Boltzmann's constant. Mo- is the coupling mass, 
M is the particle mass, and the chosen values assume a 
super-symmetric origin of the CDM particle. 

Unfortunately, the linearisation will break down at high 
redshift as the expansion-normalized scales as e~" , so the 
value by which the dark matter tilt is suppressed rela- 
tive to the baryon tilt is unclear. However, equation (|63p 
strongly suggests that the baryon-photon plasma dominates 
the momentum density at last scattering and that its tilt 
should be properly determined at 2lss by equation (|60|l with 
P = (pb + 4p.y/3)w(b) on the left-hand side. After this W(b) 
decays as 1 -I- 2 as determined by (I62|l . The usual procedure 
of assuming that all components have the same tilt under- 
estimates |u(b) I at last scattering by the ratio of the baryon- 
photon enthalpy to the total enthalpy. For the majority of 
our results, we follow the usual procedure for consistency 
with previous work. We consider the effect of the improved 
velocity analysis in Section [5.71 where we show that it will 
have a significant impact on statistical studies, but does not 
change our qualitative results. 



5.3 Parameters of the model 

For ACDM, the background model may be specified fully 
by the physical densities in CDM (ti;c,o = ^0,0/1100 with 
Ho = lOO/iioo kms~^ Mpc~^) and baryons (f^b.o) with Qa.o 
and Qk.o- The Hubble constant and 57m, are then derived 
quantities. Models with fixed luc,o and ci^b.o have the same 
early-universe history and reproduce the same acoustic peak 
structure in the CMB spectra if the angular-diameter dis- 
tance to last-scattering an d the primordial power s pectra 
are additionally held fixed (|Efstathiou fc Bond|[l999h . The 
Bianchi representation with structure constants in canonical 
form further requires us to specify e"'' , although this is of no 
physical consequence in the background. In the perturbed 
model, the current scale factor e"" determines the physi- 
cal size over which the shear eigenvectors rotate in space 
on a parallel-propagated triad. For the simplified perturbed 
model, we must additionally specify initial values for (T12 and 
ai3. Due to the rotational symmetry of the VII^ structure 
constants about ei, only m = ±1 anisotropies and polar- 
ization are generated in this model, and are proportional to 
[(a"i2T*<^i3)/^]o respectively. Varying the phase of (Ji2+jcri3 
amounts to rotating the sky about ei (reflecting the residual 



freedom in the choice of 62 and 63), while the rotationally- 
invariant content depends on (T12 -I- erf 3. We can, therefore, 
always choose ai2 = cri3 which we do for compatibility with 
previous studies. 

The morphology of the CMB anisotropy and polariza- 
tion patterns in the Bianchi model is determined largely by 
the parameters f^M.Q, ^a,o and the conformal Hubble pa- 
rameter e"°Ho. The expansion-norma lized shear (ai2/H)o 
and (ffi 3/-ff)o determine the amplitude I"*! Collins fc Hawkind 
l|l973al ) denote the conformal Hubble parameter by x, i.e. 



(65) 



where the latter relation arises directly from the FRW def- 
inition of fix.o, with K = —h in our case (Section 12. ip . 
With X, Qm,o and Qa,o fixed, variations in e°'° change phys- 
ical scales in the model (e.g. the age) but do not afi'ect the 
conformal properties. There is an approximate degeneracy 
amongst Qm.o, ^a,o a nd x that preserves the morphology o f 
the Bianchi patterns Jjaffe et al.ll2 006': 'Bri dges et al.llioOTi '). 
In our results, we follow Ijaffe et a l. (2006() by fixing Wc,o 
and ojb.o and use an ionization history consistent with these 
choices. Further specifying f2M,o and f2A,o determines Hq; 
the current scale factor is then fixed by x. 



5.4 Summary of the calculation 

We assemble a hierarchy of multipole equations for Oim, Eim 
and Bim using the results of Section U) The Thomson scat- 
tering rate f requires a model for the recombination (and po- 
tentially reionization) history, for which we use RECFAST 
(Scagcr, Sasselov fc Scott 1999). 

Starting at 2 ~ 1500, the initial power is taken to be 
zero for the polarization, with a pure dipole for the tempera- 
ture arising from the tilt of the baryon-photon fiuid. Whilst 
the universe remains optically thick, a small quadrupole 
term in the temperature and polarization distributions arises 
from the equilibrium between scattering and anisotropic red- 
shifting due to shear; this is quickly attained during the nu- 
merical integration and there is no requirement to include it 
in the initial conditions. However, we verified that starting 
significantly earlier (2 ~ 1800) made no difference to the 
final results. 

Since the power in the temperature modes declines very 
rapidly for / > 15, we truncate the hierarchy at / = 60 with- 
out any special boundary conditions. Performing the calcu- 
lation with a higher truncation {I — 120) made no difference 
to the results for / < 30. During the numerical integration 
we ensure at each timestep ^77 ^ 1, aSi] ^ 1 and tSt] ^ 1. 

We investigate the polarization properties of the 
CM B in t wo m odels on the degeneracy proposed by 
iJaffe et al.l l|2006l) . namely (a:, ^a.o, f^M.o) = (0.62,0,0.5) 
and (1.0,0.7,0.2) respectively, both with "right-handed" 
parity. The latter model is as close to a concordance value as 



The more refined treatment of tilt velocities requires one to 
specify also the fraction of baryons to dark matter and the physi- 
cal Hubble parameter. However, the Doppler terms are generally 
only a small correction to the anisotropy accumulated through 
the shear. The details of recombination introduce further depen- 
dencies on the physical densities of baryons and dark matter. 
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Figure 1. Temperature (top), i?-mode (middle) and S-modc 
polarization (bottom) maps for the Bianchi VIIjj model with 
{x, Oa.Oi f2M,o)=(0-62, 0, 0.5) and a consistent recombination his- 
tory and no reionization. The maps have been transformed to the 
observational basis {—p, eg, e^), which involves a parity change 
of the form ||4 2[I. and rota t ed to match the orientation of the tem- 
plate given in ljaffe et al.l ||2006|). The masks used in the WMAP 
polarization analvsis l lPage et al.ll2007^ are overlaid on the polar- 
ization maps. 



the Bianchi fitting allows (see Fig. 7 in lBridges et al.ll2007l ). 
In both cases, we take tJb.o = 0.022 and Wc.o = 0.110 and the 
consistent recombination history with no reionization. These 
models produce almost identical polarization patterns, for 
reasons outlined below. We briefly discuss the effects of al- 
tering the ionization history in various ways (including reion- 
ization) in Section [5.61 

In each case, we normalize such that the maximum 
temperature anisotropy corresponds to AT — ±35 /iK. 
Note that the amplitude of the polarization anisotropy sim- 
ply scales linearly with the magnitude of the temperature 
anisotropy. 




I 



Figure 2. Auto- and cross-correlation power spectra for the 
Bianchi models (x, Oa,o> ^M,o) = (0.62,0,0.5) (solid lines) and 
(x, 0: ^M,o) = (1.0,0.7,0.2) (dotted lines), normalized such 
that the maximum AT = ±35 /^K. (The units of the vertical axis 
are /^K.) The main difference between the models is a shift of 
power to larger scales (lower I) in the model with A; this is well 
understood in terms of the reduced focusing given lower f2x,o (see 
text), and causes no difference to our conclusions. Note that the 
TB correlation is negative for I < 6 and I < 5 in the respective 
models. 

5.5 Results 

The resulting temperature and polarization E- and B-mode 
maps for the flA,o = case are illustrated in Fig.[T] The level 
of the polarization is very high, approximately 1 fjK. Heuris- 
tically, this is because the shear modes considered here de- 
cay as (1 ± z)^, so that a substantial portion of the final 
temperature anisotropy can be built up between individual 
scattering events at high redshift. Because of the efficient 
conversion of _E-modes to B, (equation [52}, the B-mode con- 
tribution is of similar magnitude to the i5-mode. 
Although computing the power spectra, 

C^^ = 2i\-i o-h^cifm, (66) 

does throw away useful information in these models, it pro- 
vides a fast and efficient way to compare our results with 
known, and robust, polarization constraints. Since the mul- 
tipole hierarchy does not transfer power between different 
m values, and the implemented cosmology only generates 
anisotropics with m = ±1, in forming the power spectrum 
we are throwing away only phase information. 

Given t he position o f the Bianchi- like features on the 
sky given in Ijaffe et al.l (|2006l ). we may be confident that 
the P06, and even the P02 , mask of the WMAP polarization 
analysis (|Page et al.lbOO' ?) could not hide the polarization 
signal to a great extent (Fig. [l]). Although the relation be- 
tween the Stokes parameters and E and B is non-local, we 
note that maps of the Stokes parameters have their power 
localized in a similar way to E and B (and T). We therefore 
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Figure 3. Growth of observable r.m.s. signal in the T and E- 
and S-mode polarization plotted against a = e"""" = (1 + z)~^ 
for the models (a;, Oa o> f^M,o)=(0-62, 0, 0.5) (solid lines) and 
(x, Oa,o> f2M,o)=(l-0, 0.7, 0.2) (dotted lines). Note that the power 
grows rapidly at high redshift while the shear is still significant, 
then remains constant (although it is transferred to higher I, 
which cannot be seen in this diagram). It is for this reason that 
the polarization is remarkably strong and relatively insensitive 
to the cosmology along the line for which the VII^ temperature 
patterns are degenerate. 



calculate the full-sky power spectra without any considera- 
tion of the effect of masking nor the weighting with the in- 
verse of the (non-diagonal) pixel-pixel noise covariance ma- 
trix that were employed by the WMAP team. Given that the 
r.m.s. Bianchi signal inside the masks is lower than outside, 
we expect the effects of masking would increase the esti- 
mated Bianchi power spectra over the full-sky values plotted 
in Fig. H 

The major difference between the two parameter sets 
considered is that, for the fiA.o = 0.7 case, the distinctive 
spiral pattern is less 'focused'. This is a well understood ef- 
fect of reducing the sp atial curvature Qk.o to 0.1 from its 
original value, 0.5 (e.g. iBarrow et al.|[l985l ). and manifests 
itself as a shift of power to lower I (see Fig. [5}. The existing 
statistical studies show that distingui shing these cases ob - 
servationally is currently not possible l|Bridges et al.|[2007l ). 

There is no significant difference in the overall polariza- 
tion power. This follows because the majority of the power 
is built up rapidly at high redshifts as the universe becomes 
optically thin and the shear term has not decayed: at this 
point, the model is insensitive to the values of flA,o and 
^K,Q (Fig. [3]). Allowing u!h,o to vary introduces much more 
substantial variations in the relative level of polarization; 
however, this introduces a further degree of freedom and is 
beyond the scope of our current analysis. 

In Fig. |4l we compare the power spectra in the Bianchi 
f^A.o = model with the power expected in a 'concordance' 
model with standard, statistically-isotropic and homoge- 
neous p erturbations. The latter spectra are computed using 
CAMB (|Lewis, ChaUinor fc Lasenbvll2000l ) for two models, 




101 

/ 

Figure 4. Bianchi VII^ induced power in the CMB (solid 
lines) for (a;, Ha o, o)={0.62, 0, 0.5) and no reionization, 
compared with Gaussian power from inhomogeneities for 
(a;c,Oi '^b,Oi ""8: ^) = (0.110,0.022,0.7,0.3) with reionization op- 
tical depth T = (dash-dotted lines) and r = 0.10 (dashed 
lines). The polarization data are from the WMAP three-year 
release l|Page et al.l [20o3) . From the TE and EE power spec- 
tra alone, the Bianchi-induced polarization can mimic the ef- 
fect of early reionization in the standard scenario (the conven- 
tional interpretation of the large-scale polarization power seen by 
WMAP). However, the best-fit Bianchi model to the tempera- 
ture map clearly over-produces S-mode power compared to the 
WMAP upper limit (plotted) ruling out the simple model imme- 
diately. 



one with no reionization (dot-dashed lines) and a favoured 
reionization model (dashed lines; r = 0.10). Forming com- 
bined power spectra by adding the power from the Bianchi 
model that best fits the temperature maps to that from the 
concordance model i s inconsistent, since the models have 
different parameters (|Bridges et alJbOOTl ). However, ignor- 
ing this problem and comparing the models as 'templates' 
shows that, so far as the TE and EE power spectra are 
concerned, the Bianchi model can mimic the observed large- 
angle power that is conventionally attributed to reioniza- 
tion. Of course, the 'corrected' power in such a model would 
probably lead to an unfeasibly low estimate for r in light 
of other data such as the G unn-Peterson constraints (e.g. 
iFan. Carilh fc Keatinell2006l ). So, at least with the fiducial 
simplified dynamics outlined in Section [5. II this already pro- 
vides strong evidence against the VII/i model. 

More challenging for the Bianchi model is the B-mode 
polarization, whi ch is at a similar level to the _E-mode. In 
iPage et all (|2007l ). the _B-modes for Z < 10 are found to be 
consistent with zero with errors better than a ~ 0.1 ^K'^ 
at each multipole. At this level, the signal-to-noise on the 
B-mode spectrum in the Bianchi model should be at least 
unity for each 2 < Z < 8, and would have produced a highly 
significant detection of large-angle B-modes overall. 

Finally, the Bianchi models are not parity-invariant 
and one therefore obtains a TB and EB cross-correlation 
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(Fig. [2}. To get a rough estimate of the current statis- 
tical power of these spectra in constraining the Bianchi 
model, we compute the 'x between the model prediction 
and the WMAP estimates of Cf^ and Cf ^ available on 
the LAMBDA websit«0. We use the spectra from / = 2-16 
and, since only the diagonal errors are publically available, 
we ignore correlations between the estimates and compli- 
cations due to the shape of the low-/ likelihood. As noted 
earlier, we also ignore the effects of foreground masking and 
noise-weightinglj We find reduced values of 0.5 for TB 
and 4.3 for EB for 15 degrees of freedom. The corresponding 
figures for null ^ and Cf ^ are 0.4 and 0.6. Note that, 
although the Bianchi power is typically two orders of magni- 
tude smaller in EB than TB, the EB estimates have smaller 
errors as Cf exceeds the variance of the polarization noise 
on these scales. The interpretation of these values is that 
the data are too noisy to distinguish the Bianchi model from 
the null case for TB (both are perfectly consistent) but the 
EB spectra disfavour the Bianchi model over the null case. 

5.6 Effect of ionization history 

Since the shear decays rapidly, cr oc (1 + z)'^, our inclusion 
of a more detailed recombination calculation will affect the 
temperature maps somewhat. We take the 57a. o = model 
and run the Boltzmann hierarchy twice; first with instan- 
taneous recombination at 2; = 1100 and then with the full 
RECFAST history (with no reionization) . There is no qual- 
itative difference in the temperature maps produced, but 
there is an approximately 15 percent decrease in the temper- 
ature amplitude in the latter case. Of course, this is simply 
reflected in a slightly different estimate of {a/H)o and has no 
significant impact on previous probes of Bianchi signatures. 

However, the detailed recombination history does have 
a significant impact on the amplitude of polarization. With 
the detailed model, the amplitude is approximately five 
times larger than that derived from the instantaneous model. 
Note that this puts the amplitude of polarization in t he in- 
stanta neous model in agreement with the estimation in iReesI 
1 19681 ). It is unsurprising that the polarization is so sensi- 
tive to the recombination model, given that it arises through 
the detailed interplay of the rapidly decreasing shear and 
sharply peaked visibility function fe"'^. 

The effect of adding reionization is not as dramatic as 
for standard FRW perturbations; this is because of the high 
level of the primordial polarization relative to the tempera- 
ture signal (~ 1/25) on large scales. Adding reionization as 
early as z = 15 produces only a ~ 50 percent increase in the 
EE power. 

5.7 Effect of improved constraint model 

The numerical results presented so far have derived the 
baryon tilt assuming the same tilt for all particle species, 

^ |htt p : //lambda ■ gsf c .nasa. gov/ 1 

^ We checked that the additional variance in the power spec- 
trum estimates (when averaging over statistically-isotropic CMB 
fluctuations and noise) due to products between the two-point 
functions of the Bianchi signal and fluctuations is only a small 
correction to the errors computed by the WMAP team. 
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Figure 5. (Top) Normalized temperature map for Bianchi 
Vllft model (x, Qa,0) ^m,o)=(.0.62, 0, 0.5) with improved tilt con- 
straint; (bottom) residuals in this map compared to the standard 
constraint map (Fig. [TJ . The centre of the maps are here oriented 
down the ei axis. 




10° 10'' 



Figure 6. Comparison of power spectra with standard constraint 
(solid lines) and improved constraint (dashed lines) after renor- 
malizing the maximum temperature anisotropy to ±35 /iK. 
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consistent with previous work. In this section, we analyze the 
effect of dropping this assumption and adopting the more 
sophisticated model of Section 15.21 

One may see heuristically that the shear generally 
contributes more to the temperature anisotropy than the 
dipole, because its integrated effect (for models with mat- 
ter domination over most of the line of sight) scales as 
~ n~^o^^((j/jy")o(l + zi^ssf'^ X 0(1) whereas the dipole is 
imprinted instantaneously at scattering and scales in the 
improved model as 



VT+9h{l + ZLss) 



b,0 



HI, 



(67) 



Using this new approach with the two sets of parameters 
considered above, the ratio of the shear to Doppler contri- 
butions is ~ 6. If instead we assume the same tilt for all 
species, as in the previous section, we should replace Qh,o 
by Qm,o in equation (I67p and this ratio becomes ~ 30. 

The effects of the improved treatment of tilt are thus 
twofold: (a) an increase in the relative level of the 'distorted 
dipole' to the 'distorted quadrupole' component in the tem- 
perature maps; and (b) a decrease in the overall level of po- 
larization, given that the models are normalized to a fixed 
maximum AT/T, and only the quadrupole at high redshift 
is responsible for producing polarization. 

The revised temperature map and residual map for the 
(a;, nA,o, f^M.o) = (0.62,0,0.5) model are shown in Fig. [S] 
The power spectra are plotted in Fig. [S] As expected, the 
polarization strength is somewhat lower after renormaliza- 
tion. The effects are in accordance with our expectations: 
the difference in the temperature maps amounts to a 15% 
effect, whilst the polarization level is reduced by approxi- 
mately 10%. 

It is clear that the details of how the tilt is treated will 
impact on a detailed statistical comparison of the models 
with the WMAP data, but a full study is beyond the scope 
of the present work. However these effects are not sufficiently 
large to make the Bianchi B-mode polarization unobservable 
at the three- year WMAP sensitivity (cf. Fig. [3}, or change 
our overall conclusions. 



Il985h . although the amplitude is modified somewhat due 
to the better treatment of the ionization history. Polariza- 
tion maps, with the generality presented here, do not appear 
to have been computed before. Note also that for these, a 
detailed treatment of recombination is required for accu- 
rate results. The power spectra of our type-VIIh polariza- 
tion maps apparently put these models in cont radiction of 
the l arge-scale polarization results from WMAP (|Page et al.l 
l2007f). 



During the drafting of this paper, an analysis of uni- 
verses equivalent to Bianchi I models, tuned to account 
for the low CMB quadrupole l|Campanelli. Cea fc Tedes"col 
2006), was shown to give a sim ilar level of polarization to 
that computed here l|Ceall2007l ). This is not surprising given 
that the anisotropics are tuned to address some of the same 
problems, and that the added complications induced by the 
Vllfi geometry do not substantially alter the amplitude of 
the effect (in the simplified dynamical model) . In the type-I 
model, the temperature anisotropy and _E-mode polariza- 
tion are simply quadrupoles, and no _B-modes are produced 
(see eauation l52|l . Although the type-I model does not suffer 
the same observational constraints as type-VII;, in polariza- 
tion, the latter has the virtues in temperature of resolving 
essentially all of the large-angle anomalies. 

Our cursory glance over the available data appears to 
rule out the Vllh models employed in recent papers on the 
basis that they over-produce B-mode power. This is espe- 
cially significant given that the results hold for all models 
on the Bianchi degeneracy line given in IJaffe et al.1 l|2006l ). 
Our polarization results, combined with the failure of the 
Bianchi degeneracy region to include well-established val- 
ues for the cosmological parameters, suggest that the simple 
Vllh model, as it stands, is unsuitable to describe the avail- 
able data. However, to reject completely the hypothesis that 
our universe contains anisotropic perturbations that are ho- 
mogeneous under groups of motions with Bianchi type Yllh 
requires a fuller treatment of the dynamics of the linearized 
model (Section l5.ip . We intend to address this problem, and 
to search for statistical correlations between the morphol- 
ogy of the generalized model's polarization and the WMAP 
data, in future work. 



6 CONCLUSIONS 

We have derived the radiative-transfer equation for the 
CMB, including polarization, in all nearly-FRW Bianchi uni- 
verses in the form of a hierarchy of multipole equations 
which can be easily integrated numerically. These can be 
coupled with the dynamical (i.e. Einstein) equations to com- 
pute maps of the CMB temperature anisotropics and polar- 
ization in any such model. _B-mode polarization is generic, 
being produced in all Bianchi types except I. We applied 
these equations to the Bianchi Vllh case, with parameters 
tuned to address the anomalous features observed in th e 
CMB temperature on large scales (fjaffe et al.ll2005l , l2006h . 
Our treatment includes a more physical treatment of the 
tilt velocity in CDM models with sub-dominant baryons. 
Whilst this does not make a qualitative difference to our re- 
sults, more detailed statistical studies could well be affected 
by its ~ 20% corrections. 

Our tempera ture maps are similar to tho s e derived from 
earlier studies l|Collins fc Hawking! Il973al : iBarrow et al.l 
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